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The electronic correlations on a C20 molecule, as described by an extended Hubbard Hamiltonian 
with a nearest neighbor Coulomb interaction of strength V, are studied using quantum Monte Carlo 
and exact diagonalization methods. For electron doped C20, it is known that pair-binding arising 
from a purely electronic mechanism is absent within the standard Hubbard model (V = 0). Here 
we show that this is also the case for hole doping for < U/t < 3 and that, for both electron and 
hole doping, the effect of a non-zero V is to work against pair-binding. We also study the magnetic 
properties of the neutral molecule, and find transitions between spin singlet and triplet ground states 
for either fixed U or V values. In addition, spin, charge and pairing correlation functions on C20 
are computed. The spin-spin and charge-charge correlations are very short-range, although a weak 
enhancement in the pairing correlation is observed for a distance equal to the molecular diameter. 

PACS numbers: 71.10.Li, 02.70.Ss 



I. INTRODUCTION 

Shortly after the discovery of superconductivity in 
Ceo, it was suggested by Chakravarty, Kivelson and 
Gelfandiii 2 -^ that an electronic mechanism, in which pairs 
of electrons preferentially reside on a single molecule 
rather than on neighboring molecules, might provide the 
pairing mechanism for superconductivity. Using second 
order perturbation theory they found evidence for pair 
binding, above a threshold value of U/t rj 3. They also 
found that this attraction between doped electrons is ac- 
companied by a violation of Hund's rule, which requires 
maximal spin, for the two-electron-doped Cqq, and that 
for U/t > 3, the ground state for two-electron-doped Ceo 
has spin zero. — ^ However, recent calculations^ using 
quantum Monte Carlo (QMC) techniques, suggest that 
the repulsive Hubbard model does not lead to pairing 
on Ceo- On the other hand, there are geometries where 
pair binding is known to occur— In particular, White 
et al. in exact diagonalization (ED) studies of the ex- 
tended Hubbard model on the much smaller C12 (trun- 
cated tetrahedron) molecule have shown that a negative 
pair-binding energy (effective attraction between doped 
electrons) exists for an intermediate value of the on- 
site Coulomb interaction U [see Eq. (Q]) and Fig. 03(a)]. 
A more realistic model of the fullcrcncs would include 
longer ranged Coulomb repulsions, and it was found that 
this pairing energy also survives in C12 for modestly re- 
pulsive values of the nearest-neighbor (NN) interaction, 
V , but increasing V eventually kills the pair binding. The 
same violation of Hund's rule as in Ceo was also observed 
in C12 [see Ref. | and Fig. OH (b)]. 

With a different extended Hubbard model, Sondhi et 
al~ studied the effects of both NN interaction V and the 
off-diagonal interactions on the pair-binding energy and 
Hund's rules violation in the Cgo molecule. Using pertur- 
bative calculations, they find that the NN interaction V 
terms suppress pair binding while the off-diagonal terms 
enhance it. Goff and Phillips 7,8 considered the effects of 



both NN interaction V and longer-range terms, V, on 
the pair-binding energy, again by perturbation theory, 
and also found that the inclusion of V terms strongly 
suppresses pair binding in Ceo- 

The fact that ED studies found pair-binding for the 
smaller C12 molecule^ and the recent rapid development 
of experimental techniques for the synthesis of C20 solid 
phasesSilS, make it interesting and timely to explore cor- 
relation effects in C20, the smallest gas-phase fullerene 
molecule which has dodecahedral geometry^ In Ref. [IH, 
we briefly reported on pair-binding for electron-doped 
C20 for a wide range of values of U/t < 100, but with 
V = 0, using both QMC for U/t < 3 and ED for the 
full range of values. Using cluster perturbation the- 
or y 13 : 14 we also identified a metal-insulator transition 
near U c /t ~ 4.2 for molecular solids formed of C2o- In 
this paper, we provide further details of our numerical 
techniques and consider both electron and hole doping 
for an extended Hubbard model with both on-site and 
NN repulsion. We also study density-density, spin-spin 
and pairing correlation functions as a function of separa- 
tion on the molecule. 

The extended Hubbard Hamiltonian on a single C20 
molecule is defined as 

H = -t 2J {c\a c j<y + h - c ) + u X! n ^ ni i + V^niUj, 
(ij)a i (ij) 

(1) 

where c\ a (cj ff ) is an electron creation (annihilation) oper- 
ator on site i, indices i,j run over 20 sites of a dodecahe- 
dron, U is the on-site Coulomb interaction, V is the NN 
Coulomb interaction, and rij = n,-f is the number of 
electrons on site i. Our goal is here to focus on strong cor- 
relation effects in C20 using exact numerical techniques. 
The Hamiltonian Eq. ([T]) is a simplified model of C20 
but it still largely captures such correlation effects. We 
calculate ground state energies as a function of both U 
and V for neutral, one- and two-electron dopings. Com- 
parisons among these energies show that the electronic 
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FIG. 1: Huckel molecular orbitals of a neutral dodecahedral 
C20 molecule. — 



pair-binding energy A b (21) = £(20) + £(22) - 2£(21) 
is positive (repulsive) for the parameter ranges studied 
(0 < U/t < 3 for V/t = 0.2 and 0.20 < V/t < 0.46 for 
U/t = 1). This implies that it is energetically favorable 
for two electrons to stay on different molecules as opposed 
to the same molecule. We also find that the existence of a 
NN Coulomb interaction V enhances this tendency, as ex- 
pected, in order to reduce the intramolecular Coulomb in- 
teraction energy. For hole doping, the corresponding hole 
pair-binding energy A b (19) = £(18) + £(20) - 2£(19) 
is again positive (repulsive) for the parameter range 
(0 < U/t < 3 and 1^ = 0), i.e., there is an effective repul- 
sion between two doped holes on the same C20 molecule. 

Unlike the case of Cgo, the highest occupied molecu- 
lar orbital (HOMO) of the neutral C20 molecules, in the 
weakly interacting limit, is a four-fold orbitally degener- 
ate level occupied by two electrons. (See Fig. [TJ) Hund's 
rules predict for this case that the two electrons occupy 
different orbitals and have total S = 1, implying that, 
in the absence of a Jahn- Teller distortion, the neutral 
molecule has a magnetic moment. In previous worki^, 
for V = 0, we have confirmed this magnetic moment for 
< U/t < 3 and shown that at the metal-insulator tran- 
sition, U c , the ground-state changes from a spin triplet 
to a singlet for neutral C20 and from S = 2, through 
S = 1, to S = for Clo . Here we extend this analysis 
to determine ground state spin configuration for neutral 
C20 for a fixed value of U/t = 2 as a function of V, and 
find a level crossing between V/t = 1 and V/t = 1.5 for 
spin triplet and singlet states. For U/t = 2, we estimate 
the critical V c /t to be 1.1 for the spin triplet to singlet 
transition of the neutral molecule. In light of our results 
for V — 0, we expect that, in this case too, the mag- 
netic transition at V c will coincide with a metal-insulator 
transition for molecular solids formed of C2o- We also in- 
vestigate the pair-binding energy for the hole doped case 
for both V = and V 7^ 0, and examine the effect of a 
non-zero V on Hund's rule. 

The occurrence of orbital degeneracy and the resulting 
magnetic moment are tied to the icosahedral symmetry 



of the molecule. Simple molecular orbital calculations 
strongly suggest that the molecular symmetry is lowered 
by a Jahn- Teller effect from Ih to £3^, with the HOMO 
being a non-degenerate singlet. ~ However, the correla- 
tion effects that give rise to Hund's rule compete with this 
tendency to form a singlet ground state, and hence they 
also compete with the Jahn- Teller effect. As reported 
previously^ we find that when the on-site Coulomb in- 
teraction U/t is sufficiently large (U > 4.2t), the ground 
state is gapped with S = and the Ih symmetry is likely 
stable against a £3^ distortion. In order to more exclu- 
sively focus on the effects of the non-zero V term we shall 
here assume that the icosahedral symmetry is unbroken 
even for smaller U values. 

In the next section, we briefly introduce the projection 
quantum Monte Carlo (PQMC)^i and ED methods for 
this model. This is followed, in Section ITTTl by a compar- 
ison of PQMC with ED results on a C12 and a discussion 
of Hund's rule violation in C12. In section HVl we focus on 
the C20 molecule. Hole pai-rbinding in C20 is discussed 
and the influence of a non-zero next nearest neighbor V 
on the pair-binding is investigated and results for the 
triplet-singlet transition with V are described along with 
calculations of several correlation functions in the C20 
molecule. Section IVl contains discussion and conclusions. 



II. METHOD 
A. PQMC 

As noted in Ref. [IH, the idea in PQMC simulations of 
the extended Hubbard model is to decouple the two-body 
interaction terms (both U and V terms) in the parti- 
tion function by means of discrete Hubbard-Stratonivich 
transformations. — The resultant one-body terms are 
coupled to several auxiliary Ising spin fields that live ei- 
ther on the lattice sites (U term) or on the lattice bonds 
(V term). One such discrete transformation in the V 
term is given by 

-ATVni„n tf _ }_ry \2<?°'f(ni c ,-n jl 3)-^^(ni a +n j/3 ) 

e ~ 2 Kf} e 

(2) 

where a,/3 =|,i, crff = ±1 is the auxiliary Ising 
spin on bond («,j), At is the discrete imaginary time 
slice in PQMC, and the parameter A2 is determined by 
tanh 2 (A2/2) = tanh(^j^). The same decoupling equa- 
tion applies for on-site Coulomb interactions, i.e., i = j, 
except that the constant V is replaced by U and A2 by 
Ai, which is similarly given by tanh 2 (Ai/2) = tanh(Al^i). 
These one-body fermionic terms in the partition function 
can then be explicitly traced out, leaving traces over the 
auxiliary Ising spins, which can be evaluated by Monte 
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Carlo (MC)±2. 

Z = J2H^ + BLi^B^ia) ■ ■ ■ Bt(c 

{a} a 

= ^detO({<r}) T detO({cT});, 



(3) 



where {cr} = {cr 1 , cr 2 , cr 3 , a 4 , a 5 } is the set of five species 
of Ising fields, with a 1 representing the on-site Ising spins 
and cr 2-5 the NN bond Ising spins (one for each of the 4 
spin configurations). The Bi matrices are defined as 



((a) = e -ArJf/3 e l^(0 e -ArJf/2 j 



-i for i,j NN, 
otherwise, 



(4) 
(5) 



W*(l) = a[* y Ai*< (0 + 6 {ij) \ 2 ]T <], (6) 



1 for ij NN, 
otherwise, 



(7) 



where £ = 1, • • • , L is the time slice index, and a = ±1 
denotes the two determinants in Eq. (|3|). 

A complete MC sweep through the lattice will there- 
fore consist of trial flipping of one species of auxiliary 
Ising spins on all the lattice sites and trial flipping of four 
species of auxiliary Ising spins on all the NN bonds in the 
lattice system. Fast calculation of the probability ratio 
in flipping one bond Ising spin at one time slice is still 
possible using the local update technique, 20 except that 
one needs to apply the probability ratio formula twice for 
each bond Ising spin flip (which affects two sites). 

We remark that in this decomposition scheme it is 
possible to treat even longer range Coulomb interactions 
[e.g., next nearest neighbor (NNN) Coulomb interactions, 
etc.] by introducing more species of auxiliary Ising spins 
that live on these longer bonds. The only problem is 
that one needs to walk through a larger and larger phase 
space of the auxiliary Ising spins during the MC sim- 
ulations, which will, of course, increase the computa- 
tion time. Practically, we find that, to collect the same 
amount of data, the CPU time doubles for V ^ com- 
pared with the V — case. 

In a typical calculation the projection factor (3 in 
PQMC was taken to be j3 = 10/t, and the discrete time 
slice was set at At = 0.05 /t. 10 3 MC warm-up sweeps 
through the whole space-time lattice are typically per- 
formed before collecting data. To estimate the statistical 
errors, we use the same method as was used in Ref. |4|. 



B. Exact Diagonalizations (ED) 

The exact diagonalizations on C12 are done using stan- 
dard Lanczos techniques and we therefore focus on the 
ED of C2o- We always use total particle number N and 
total S z as quantum numbers since they are conserved 




FIG. 2: (Color online) Dodecahedral C20 geometry in 2D 
view. Solid and empty points denote two sets (orbits) of car- 
bon atoms divided by the S10 symmetry. 



and we perform ED in the corresponding reduced Hilbcrt 
space. In addition, ED are performed using the Sw sub- 
group symmetry present in the point group J/j. The im- 
proper rotations generated by the elements of <Sio can 
be visualized as a rotation of an angle 27r/10 around the 
center of a pentagon followed by a reflection in a plane 
perpendicular to the rotation axis. This is illustrated 
in Fig. [2] where the numbering of the sites is to be un- 
derstood in the following way: the sites 1 through 10 
are shifted up by 1 (modulo 10) under Sw and the sites 
11 through 20 are shifted in a similar manner. Hence, 
under the action of the Sw group two different orbits ex- 
ist, marked by the solid and open points in Fig. O Many 
other symmetries exist but the Sw symmetry is large and 
relatively easy to implement, and we have not exploited 
additional symmetries since the added cpu-time needed 
to implement them was significant enough to offset the 
time gained from reducing the size of the Hilbert space. 
The Sio quantum number can be thought of as a pseudo 
angular momentum, jw, and for each value of N and S z 
we have to find the value of jw that corresponds to the 
ground-state. In many cases it is not an obvious value 
and it is often non-zero. In the accompanying tables we 
show the values of jio corresponding to the listed ener- 
gies and in Table IIVI we show complete dispersion of the 
lowest magnetic modes for neutral C20 as a function of 
jio- 

The calculations are fully parallelized Lanczos calcu- 
lations executed on SHARCNET computers. A typical 
calculation performed at half-filling for N — 20, Sz = 
that, after Siq symmetry reductions, requires a Hilbcrt 
space of M = 3, 418, 725, 024, is performed with P = 64 
cpu's using about 540 seconds of cpu-time (for each cpu) 



4 



per Lanczos iteration. The memory requirement for this 
example is roughly 2.1Gb per cpu. Excellent convergence 
is always observed with less than 300 Lanczos iterations, 
typically less than 200. 

The heart of the Lanczos calculation is the matrix vec- 
tor multiplication that in this case has to be implemented 
in parallel. As one of several choices, we have chosen to 
have each cpu apply the full matrix to one section of the 
vector with each cpu returning the corresponding section 
of the resulting vector. The partial results from each cpu 
therefore needs to be communicated between all P pro- 
cessors with each processor communicating to all others. 
Due to the size of the involved Lanczos vectors ( 40-60Gb) 
which greatly exceeds the available per-cpu memory, it 
is necessary to repeat this P x P communication step 
many thousands of times per Lanczos step. The commu- 
nication step therefore quickly becomes the bottle-neck 
in the calculation unless it can be done very efficiently. 
Fortunately, this is possible using non-blocking commu- 
nications where the individual cpu's do not wait for a 
communication to complete. The draw back of using non- 
blocking communications is that buffer space has to be 
allocated until it has explicitly been verified that the com- 
munication has been completed. We have implemented a 
dual buffer strategy yielding an extremely efficient com- 
munication step. The cpu-time spent per cpu is, for all 
accessible number of processors we have been able to 
check, overwhelmingly dominated by actual calculations 
rather than communications. For a fixed TV we have then 
observed almost linear scaling for P = 64,128,256,384 
and 512. The great advantage of this approach is that 
the complexity of Lanczos calculations scale with the size 
of the HUbert space, TV, as TV log TV. Neglecting the loga- 
rithm, a doubling of the size of the Hilbert space, TV, can 
then be almost compensated by doubling the number of 
processors P. 
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-9.4647669965 -9.466(2) 0.97 
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1/2 1/2 -6.8287003500 -6.829(4) 0.33 


E 13 


3/2 3/2 -6.0844214907 -6.059(6) 0.20 


Ei4 


-4.1568425864 -4.11(1) 0.11 


Ei4 


1 1 -4.0772924523 -4.080(5) 0.34 


Ai.o 
Ai.o 
A b (13) 


(1/2,0) 2.6360666465 2.637(4) 
(3/2,0) 3.3803455058 3.407(6) 
(0,0,1/2) 0.0357911171 0.08(1) 



TABLE I: Comparison of ED and PQMC calculations on the 
truncated tetrahedron (12 sites) at U = 2t and V = 0.2t. 
E n (Sz) is the energy of a system with n electrons and z- 
component of total spin S z . A„ im is the energy difference 
E 12+n (S?) - E 12+m (ST) with (S",ST) given in the second 
column. For binding energies Ab(n) the second column shows 
(Sz +1 , S™- 1 , S") - the S z values for 3 states involved in its 
calculation^ 
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FIG. 3: (Color online) (a) Variation of the pair-binding energy 
A b (13) = £(12) + £(14) - 2.25(13) of a truncated tetrahedron 
molecule (C12) with U and V as in Fig. 3 in Ref.i. (b)Hund's 
rules violation in the two-electron doped C12 molecule, where 
A£(14) = ^(triplet) - Eia (singlet). 



III. RESULTS FOR C 12 

Before turning our attention to the C20 molecule we 
investigate the simpler C12 molecule in the truncated 
tetrahedron configuration. As mentioned above, previ- 
ous studies^ have found a negative pair-binding energy 
on this molecule that, however, became positive (repul- 
sive interaction) in the presence of a sufficiently large 
V. The purpose of this investigation is two-fold. First 
of all, we want to verify the correctness of our numeri- 
cal approach while at the same time highlighting some 
of the subtleties of interpreting the PQMC data. Sec- 
ondly, due to the relative ease with which calculations 
can be performed on this molecule it allows for a rather 
detailed study of the correlation between the negative 
pair-binding energy and a violation of Hund's rule for 
the two electron doped molecule^. 



A. Tests on the C12 molecule 

To test our ED program, we use the same parameters 
as in Ref. [5j and we are able to reproduce the same pair- 
binding energy as shown in Fig. [3] (a). In Table HI we see 
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good agreement between PQMC and ED energy values 
within statistical error bounds. An exception is found for 
E14 and S z = 0, where the PQMC result lies a bit higher 
than the ED energy value. This is due to the mixture of 
singlet and triplet components in the S z — sector and 
the near degeneracy of these two states that makes the 
projection of the singlet ground state out of the mixed 
state difficult X We will see that this difficulty does not 
occur for C20, where the ground state with two-electron 
doping is in the spin-2 sector for U/t < 3. Hence the pair- 
binding energy extracted for C20 by PQMC for U/t < 3 
is more reliable than the one for C12. 

B. Hund's rule violation for C^J 

In the perturbation theory studies of pair-binding in 
the larger fullcrcnc Cgo, 1,2 it was noted that a negative 
pair-binding energy (effective attraction) was correlated 
with a violation of Hund's rule for the two-electron doped 
molecule; i.e., that for Cg^~ , the ground state was found to 
be a singlet. Although our QMC results did not support 
the existence of pair-binding in Cg^~ and found a spin- 
triplet ground state, it is of interest to examine the corre- 
lation between pair-binding and the violation of Hund's 
rule in C12. The non-interacting V = U = neutral 
molecule has completely filled levels and hence a total 
spin zero. Added electrons therefore enter an unfilled 
level with an orbital degeneracy of 3. Hund's rule would 
then predict C^ to have total 5 = 1. What we find 
is that the ground state of C 2 ^ is a singlet both when 
the pairing is attractive and when it is driven repulsive 
by increasing the nearest neighbor repulsion V. This is 
shown in Fig. [3fb) where the singlet state is found to 
lie below the triplet state for both positive and negative 
pair-binding energies, for the range of U and V studied. 
Thus, for this case, Hund's rule is found to be violated 
where pair-binding occurs as well as where it does not. 



IV. RESULTS FOR C 20 MOLECULE 

We now turn to the more interesting case of the C20 
molecule. Compared with the V = case, where PQMC 
already has a sign problem for the non-bipartite dodeca- 
hedral molecular geometry, the NN Coulomb interaction 
V terms introduce more sources of negative probability 
weight, lowering the average value of the sign. Fig. [4] 
shows the average sign for these two cases. For the worst 
case (N = 20, S z = l,U/t = 3, V/t = 0.2), where the 
average sign is as low as 0.05, we have collected 7.2 x 10 7 
MC lattice sweeps. This gives a relatively large but nev- 
ertheless meaningful error bar. [See Fig. O (a).] For 
other parameter values, we have collected about 2.2 x 10 7 
MC sweeps. The acceptance ratio for the on-site Ising 
spin trial flipping ranges from 80% (U/t — 3) to 93% 
(U/t = 1), while that for the bond Ising spins is about 
95% due to the small value of V/t = 0.2. 
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FIG. 4: Average sign behavior for both V/t — (solid sym- 
bols) and V/t — 0.2 (hollow symbols) at different fillings 
N = 20,21,22. The lines connecting the points are guides 
to the eye only. 



A. Pair-binding energy 

Table Hi] shows the energies of the C20 molecule at 
different fillings from PQMC and ED for U/t = 2 and 
V = 0. Both ED and PQMC predict the ground states 
to be in the same spin sectors for the molecule, and the 
calculated energies are in agreement within MC error 
bounds. 

In order to understand the comparison of PQMC and 
ED data in TableHIl it is important to recognize a system- 
atic weakness of PQMC which is that, when the ground 
state is a spin multiplet, the different partners appear 
to have different energies, increasing with decreasing val- 
ues of \S Z \, because the states with smaller values of l^j 
mix with higher lying states that have the same value of 
liS^]. In general except for statistical error, a state with 
S z = will appear to lie above its partners with the same 
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TABLE II: Comparison of ground state energies from ED and 
PQMC calculations on the C20 molecule at U/t = 2 and V = 
0. See the caption in Table U for the corresponding definition 
of various quantities. 
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FIG. 5: (Color online) Electronic pair-binding energies 
A b (21)/t as a function of U/t and V/t from ED and PQMC 
simulations, (a) The variation of pair-binding energy with 
U/t for fixed V/t values, (b) The variation of pair-binding 
energy with V/t for fixed U/t = 1. The lines connecting MC 
and ED points are guides to the eye only. 



total 5. This tendency is apparent in the results for E20 
(5=1), E 2 i (5=3/2), and E 22 (5=2). Conversely, if a 
ground state with S z = lies below a state with S z = 1, 
we expect the ground state to be a singlet. However, in 
this case, the value of the ground state energy will be 
perturbed upward by any admixture of the next higher 
state with S = 1, S z = 0, as happened for Eu(S = 0) in 
Table Q] In general it is also true that accurate PQMC 
results are more easily obtained when the average sign is 
close to 1 compared to when the average sign is small. 

Pair-binding energies Ab(21)/i (electron) and 
Ab(19)/t (hole) as a function of both U/t and V/t are 
shown in Fig. [5] and [SI respectively. For V = 0, for both 
electron and hole doping, we see that the pair-binding 
energy is always positive (repulsive) for U/t > 0, and 
increases with increasing U/t. This is the same behavior 
as we observed for the Ceo molecule X Turning on the NN 
Coulomb interaction V {V/t — 0.2 in Fig. [5^a)) increases 
the pair-binding energy further. Hence, putting two 
extra electrons on the same neutral molecule becomes 
more costly when the NN Coulomb interaction is not 
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FIG. 6: (Color online) Hole pair-binding energies Ab(19)/t as 
a function of U/t for V — from ED and PQMC simulations. 
The lines connecting MC points are guides to the eye only. 



negligible. Panel (b) in Fig. [5] shows the variation of 
the pair-binding energy as a function of V/t for fixed 
U/t = 1. Again the pair-binding energy is positive 
(repulsive), and generally increases with V/t. The 
agreement between ED and PQMC results is fairly good 
and even though the PQMC data show some tendency 
to non-monotonic behavior for this interval of V/t, the 
ED results show that this is explained by the natural 
statistical spread of the data. Hence, in the regime 
V < V c , U <U C , the pair-binding energy increases with 
both U and V and energetically, it becomes increasingly 
favorable for two electrons to stay on two different C20 
molecules. However, we note that, for V = 0, U > U c 
it was previously founcU^ that the pair-binding energy 
decreases with U, reaching a minimum at U/t ~ 10, 
before increasing and reaching a finite value in the 
U — » 00 limit. 





S S z 


U/t = 3 


S S z 


U/t = 5(ED) 


jio 


E20 


1 


-17.04(2) 





-12.111284292 


5 


E20 


1 1 


-17.036(6) 


1 1 


-11.877033283 


0, ±2 


E21 


3/2 1/2 


-15.29(6) 


1/2 1/2 


-9.1165560273 ±1,±3,5 


E21 


3/2 3/2 


-15.529(5) 


3/2 3/2 


-8.9633623599 


±1,±3 


E22 


2 


-13.936353 


1 


-5.9715313615 


±2, ±4 


E22 


2 1 


-13.81(2) 


1 1 


-5.9715313615 


±2, ±4 


E22 


2 2 


-13.935(1) 









TABLE III: Ground state energies for neutral, one- and two- 
electron-doped C20 molecules at U/t = 3,5 and V = from 
PQMC and ED, which shows a transition between Hund's and 
anti-Hund's states at 3 < U/t < 5 for neutral, one-, and two- 
electron-doped molecules, respectively. Data without error 
bars are from ED. 



7 



u s 


jio = jio = ±1 jio = ±2 jio = ±3 jio = ±4 jio = 5 


2 
1 1 


-20.5920234655 > -19.90 -20.5920234655 > -19.90 -20.5920234655 -20.0527029539 
-20.5983834340 -19.9776970001 -20.5983834340 -19.9776970001 -20.5981592741 -19.9634427213 




-12.0123014488 -11.6726562451 -12.0123014488 -11.6726562451 -12.0123014488 -12.1112842959 
-11.8770332831 -11.8472120431 -11.8770332831 -11.8472120431 -11.8103044760 -11.8118179567 




-8.0452584717 -7.806831 -8.0452584717 -7.8068365859 -8.0452584717 -8.1803385740 
-7.9497836200 -7.9415479844 -7.9497836200 -7.9415479844 -7.8490047592 -7.9156714009 



TABLE IV: ED results for the dispersion of the lowest singlet and triplet states with jio for neutral C20 with U/t — 2, 5 and 8 
in all cases with V/t = 0. 



B. Hund's rule 

It is also clear, from the data in Table [TT1 and IIIIl that 
Hund's rule is obeyed for the corresponding range of pa- 
rameters, i.e., U/t < 3, V = 0. That is, the ground 
state for 20 through 22 all have the maximum values 
of total spin for electrons outside the C%q core, ranging 
from total spin 1 for 20 electrons through total spin 2 
for 22 electrons. This behavior occurs in the range of 
parameters where PQMC converges (for maximal |<S 2 | as 
discussed above.) As U/t is increased above 3, the sign 
problem prevents reliable PQMC calculations. This diffi- 
culty does not arise in ED where accurate calculations are 
possible for essentially any value of U/t. We have used 
ED to explore what happens for larger values of U/tJ2. 
For example, results for U — 5t are shown in the right 
hand columns of Table IIIIl Here Hund's rule is clearly 
violated. For 20 electrons, the ground state has spin zero; 
for 21 electrons the ground state has spin 1/2; while for 
22 electrons the ground state has spin 1. Clearly there 
are level crossings in the range 3 < U/t < 5. Additional 
results in this regime are given in Ref. [H]. ED also al- 
lows the calculation of the spin gap, the gap between the 
ground state and the lowest lying excited state with dif- 
ferent total spin. Results are shown in Table IIVI for a 
neutral C20 molecule with U = 2,5,8 and V — 0. When 
the metal-insulator transition occurs in the vicinity of 
U c /t ~ 4.2, the ground-state spin changes from an or- 
bitally degenerate S — 1 for U < U c to a non-degenerate 
singlet for U > U c . From the results presented in table llVI 
we see that it is the singlet state at jio = 5 that moves 
toward the bottom of the spectrum with increasing U 
and eventually, for U > U c becomes the ground-state. 
Focusing on the case V/t = 0, we see from table IIVI that 
at U/t — 2, the ground-state energy is a singlet E 1 jt = 
—20.5983834340 with a gap to the lowest lying singlet of 
AE 1 '°/t = 0.0063599685. Here the superscripts denote 
the spin of the ground- and excited states, respectively. 
For U/t > 5 wc find that the ground-state for the It con- 
figuration now is a non-degenerate singlet, S — 0, with 
energy E°/t = -12.1112842922. The lowest lying triplet 
excitation with AE°' 1 /t = 0.2342510092. This picture 
continues to hold for larger U/t with the triplet gap at 
U/t = 8 only slightly larger, AE°' 1 /t = 0.2305549540. 

Next we explore the ground state spin of the neutral 
molecule with different V/t values for a fixed U/t = 2. 



Using ED techniques we determine that the ground-state 
for V/t — 1 and V/t — 1.5 in both cases occur for 
jio = 0. However, the ground state changes from a 
spin triplet for V/t = 1 to a spin singlet for V/t = 1.5. 
Specifically, we find at V/t = 1, E(Singlet)=5.702018 
and E(triplet)=5. 639496, whereas for V/t — 1.5 we find 
E(Singlet)=17.318536 and E (triplet) =17.499741. By as- 
suming a linear dependence of the energy on V/t in this 
region, we determine that the level crossing occurs near 
V c /t ~ 1.1 for U/t = 2. 



C. Correlation functions 

We have also investigated what other correlations 
might be induced in the C20 molecule by calculating the 
following correlation functions: charge-charge, spin-spin, 
and pairing correlations as a function of lattice distance. 
Similar calculations for the Cqo molecule have been re- 
ported in Refs. [2l] andl22l. 

We define the correlation functions with respect to lat- 
tice site 1 in the neutral molecule: (nirii) is the charge- 
charge correlation, (Si ■ Si) is the spin-spin correlation, 
and (c\ a c\ _ .Cj j _ cr Ci )Cr ) is the pairing correlation, where 
i = 1, . . . , 20. In Fig. [7J we show the variation of these 
correlations for U = 2t and 3t as a function of lattice 
spacing du/R, where du is the distance between site 1 
and i, and R is the molecular diameter. For the dodeca- 
hedral geometry, there are only 5 inequivalent neighbors, 
all at distinct distances. One can understand the on-site 
correlations in terms of the probabilities, p n , n=0,l,2, for 
having n electrons on each site. Then, the on-site corre- 
lations functions are (n\) = pi + 4^2, (Sf) = 3pi/4, and 

( C L C l- CT C l-<7Cl, ff ) =P2- 

In Fig. UJa) we show results for the charge-charge cor- 
relation function for 2 different values of U with V = 0. 
As expected, the on-site charge-charge correlation is re- 
duced by an increase of the on-site Coulomb interaction 
U [panel (a)]. At larger distances, the charge on site 1 
and i are uncorrelated. The unit value of the charge- 
charge correlation corresponds to uniform distribution of 
charge. 

Fig.[7jjb) shows the spin-spin correlation function again 
for U/t = 2,3 with V = 0. The NN spin-spin correlation 
has a negative finite value, and its magnitude is enhanced 
by a larger U value. For spatial distances larger than 1 
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FIG. 7: (Color online) Variation of (a) charge-charge, (b) 
spin-spin, and (c) pairing correlation functions for U = 2t, 3t 
and V — for a C20 molecule with respect to the lattice site 
spacing, da is distance between site 1 and i. R is diameter 
of C20 molecule. 



FIG. 8: (Color online) Variation of (a) charge-charge, (b) 
spin-spin, and (c) pairing correlation functions for U = 3t 
and V = 0,0. It, 0.2t for a C20 molecule with respect to the 
lattice site spacing. 



we see that this correlation function quickly approaches 
0. Similar behavior has been observed for the spin-spin 
correlation in the Ceo molecule in Ref. HH and it was sug- 
gested that the rapid decay of the spin-spin correlation 
function was indicative of a resonant valence bond (RVB) 
or "spin dimer" state. The similarity between our results 
and those of Ref. Hi], suggest that the spin correlations 
in the ground-state of C20 also might be described by 
considering valence bond states including only dimers of 
relatively short length. 

QMC results for the pair correlation are shown in 
Fig. Efc) with U/t = 2, 3 and V = 0. Interestingly 
there is a peak of the pairing order when site 1 and i are 
NN sites. This again supports the RVB or "spin dimer" 
model for the ground-state. Beyond the nearest-neighbor 
distance, the pairing correlation function, along with the 
other correlation functions, is very close to its uncorre- 
cted value except for du = R, where the pairing order 
parameter is slightly enhanced. At the same time the 
spin-spin correlation is slightly negative showing an anti- 
ferromagnet correlation. This corresponds to the "dumb 
bell" model proposed in Ref. [22J where electron pairs are 



formed at the maximal distances of the molecular diam- 
eter. We note that in the present case, the enhancement 
of the correlations at the distances of R corresponding to 
this "dumb bell" pairing is relatively weak. 

We have also studied the influence of a non-zero V on 
the correlations. In Fig. [5] we show results for a fixed 
U/t = 3 and three different values of V = Q, O.li, 0.2i. 
Clearly, the effect of the NN Coulomb interaction V on 
these correlation functions is relatively weak, with the 
curves being almost identical for the range of V consid- 
ered here. 



V. CONCLUSIONS 

In this paper we have studied the extended Hubbard 
model on a C20 molecule through ED and PQMC simu- 
lations. The comparison clearly elucidates the relative 
strengths of the two methods. PQMC is possible for 
much larger systems than can be treated by ED. How- 
ever, ED has been applied successfully to the Hubbard 
model on 20 sites with 18-22 electrons, by making effec- 
tive use of the capabilities of a large number of coupled 
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processors. PQMC works best when the ground state is 
well separated from excited states with the same value 
of S z . As a result, ground states with larger total spin S 
and maximal \S Z \ are most accurately determined, while 
ground states with S = are sometimes problematic. 
This behavior was also found in our earlier work on C^o j* 
and the comparison of ED and PQMC results for C20 is 
consistent with and lends confidence to those earlier re- 
sults. 

The pair-binding energy for C20 shows that ex- 
tra added electrons (holes) prefer to sit on different 
molecules, rather than to reside in pairs on molecules. 
This rules out the possibility that the extended Hub- 
bard model on a single C20 molecule can produce an ef- 
fective attraction between electrons (holes) from purely 
electronic interactions. Our earlier work showed that this 
conclusion applies to the Ceo molecule as well.— We also 
find that Hund's rule is obeyed for U/t < 3 and small 
values of V and that larger values of U and V lead to 
level crossings and ground states for which Hund's rule 
is violated. For fixed V = 0, we have determined that 
this transition happens between U/t = 3 and U/t = 5, 
at U c /t ~ 4.2. And for fixed U/t — 2, as a function of V, 
we have determined that this transition happens between 
V/t = 1 and V/t = 1.5, at V c ~ 1.1. As was the case at 
the transition occurring at U c /t ~ 4.2 for V = 0, we 



expect this transition to coincide with a metal-insulator 
transition for molecular solids formed of C2o- More gen- 
erally, for U/t < 3 and V/t < 0.2, we find that the spin, 
charge and pairing correlations fall off rapidly even in the 
presence of NN Coulomb repulsion. It is an interesting 
open question if molecular solids formed of C20, in par- 
ticular away from half-filling, would display non-trivial 
order for V > V c . The answer to this question would be 
numerically demanding and we have therefore left it for 
future work. 
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